Relativistic Compact Objects in Isotropic Coordinates 

M. K. Malf and T. HarkcQ 

Department of Physics, The University of Hong Kong, 
Pok Fu Lam Road, Hong Kong SAR, P. R. China 
(Dated: February 7, 2008) 

Abstract 

We present a matrix method for obtaining new classes of exact solutions for Einstein's equations 
representing static perfect fluid spheres. By means of a matrix transformation, we reduce Einstein's 
equations to two independent Riccati type differential equations for which three classes of solutions 
are obtained. One class of the solutions corresponding to the linear barotropic type fluid with an 
equation of state p = 7p is discussed in detail. 
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I. INTRODUCTION 



Relativistic stellar models have been studied ever since the first solution of Einstein's field 
equation was obtained by Schwarzschild for the interior of a compact object in hydrostatic 
equilibrium. The search for the exact solutions is continuously of an interest to physicists. 
The most general exact solution for Vaidya-Tikekar isentropic superdense star was obtained 
by Gupta and Jasmin ]| . Hernandez and Nunez presented a general method for obtain- 
ing static anisotropic spherically symmetric solutions satisfying a nonlocal equation of state 
from known density profiles. Several new exact solutions for anisotropic stars of constant 
density were presented in Dev and Gleiser S. In Q an algorithm was proposed to gen- 
erate any number of physically realistic pressure and density profiles for spherical perfect 
fluid distributions without evaluating integrals. The gravitational field equations for static 
stellar models with a linear barotropic equation of state and with a polytropic equation of 
state p = kp l+l l n were recast respectively into two complementary 3-dimensional regular 
system of ordinary differential equations on compact state space 0,0]. These systems were 
analyzed numerically and qualitatively, using the theory of dynamical systems. Schmidt 
and Homann [t| discussed numerical solutions of Einstein's field equation describing static 
spherically symmetric conglomerations of a photon star with an equation of state p = 3p. 
Upper limits for the mass-radius ratio were derived for compact general relativistic objects 
in the presence of a cosmological constant in Q| and for the charged general relativistic fluid 



spheres in The effect of a cosmological constant on the structure of the general rela- 
tivistic objects was analyzed by Bohmer [ljj]. The equations describing the adiabatic, small 
radial oscillations of general relativistic stars have been generalized to include the effects 
of a cosmological constant in jllj]. Interior solutions for charged and neutral anisotropic 



fluid spheres, satisfying all the required physical conditions, were obtained in |12j and 13 1. 
Excellent reviews of the relativistic static fluid spheres were given in Finch and Skea 1J| , 
Herrera and Santos j3| and Delgaty and Lake jig . 

For a static fluid configuration Einstein's equations represent an under-determined system 
of nonlinear ordinary differential equations of the second order. For the special case of a 
static isotropic perfect fluid, the field equations of Einstein's theory can be reduced to a 
set of three coupled ordinary differential equations in four unknowns. To obtain a realistic 
stellar model, one can start with an equation of state. Such input equations of state do 
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not normally allow for closed form solutions. In arriving to exact solutions, one can solve 
the field equations by making an ad hoc assumption for one of the metric functions or the 
energy density, hence the equation of state being computed from the resulting metric. 

It is the purpose of the present paper to present a matrix method for obtaining new 
classes of exact solutions for Einstein's equations representing static perfect fluid spheres. 
By means of a matrix transformation, we reduce Einstein's equations to two independent 
Riccati's differential equations for which three classes of solutions are obtained. One class 
of the solutions corresponding to the linear barotropic type fluid with an equation of state 
p = 7p is discussed in detail. 

The present paper is organized as follows. In Section II we reduce the basic equation 
describing the interior stellar evolution to a matrix equation which is separated into two 
independent Riccati's differential equations. Three classes of solutions of the gravitational 
field equations are presented and the stiff solution is discussed in Section III. We conclude 
our results in Section IV. 



II. FIELD EQUATIONS AND THE MATRIX METHOD 



For a relativistic fluid sphere model the line element inside the matter distribution takes 
the form 

ds 2 = e^dt 2 - e A « (dr 2 + rW) , (1) 

where dVL 2 = d9 2 + sin 2 9d(p 2 . 

Assuming that the matter content inside the spheres is a perfect fluid, the field equations 



take the form 
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By introducing the transformation [1 



R = In r, 



(5) 



equations ©-© become 
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where the dot denotes The set of coordinates resulting after the re-scaling of the radial 
coordinate r by means of the transformation (JHJ) are called isotropic coordinates. 

From the isotropic pressure condition, we obtain the following basic equation describing 
the structure of general relativistic stellar type objects in isotropic coordinates: 
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By means of the following transformations 



- Az> - 2 [v + A 



0. 



(9) 



A 



udR, v 



wdR, 



(10) 



and with the use of the two matrices 

L = [(2u - 4m) + {2w - Aw)} , 
Equation can be written as a matrix equation in the form 

L = M T AM. 



(11) 
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We introduce now two new variables -u_ and u + which are obtained by means of a linear 
transformation, described by the matrix K and applied to the matrix N, 
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In the new variables, equation ()12|) reduces to the form L = N T K T AKN . We choose 
the elements of the matrix K so that 



K T AK 



-y/2 
y/2 

As a result of these transformations the right hand side of the matrix equation (|T2|) can be 
diagonalized, and the resulting equations satisfied by the new unknown functions u_, u + are 

(111— q dll.i Q t-^ / / , \ 

ra_— — — 2m_M_ + u_ = —m + —— + 2m + u + + u + = F (R) , (14) 
dR dR 



where m T = y 2 =p y/2 and a solution generating function F (R) is introduced. 



III. THREE CLASSES OF SOLUTIONS OF THE STATIC FIELD EQUATIONS 

In the previous Section, by means of a matrix transformation, the cross-term \v cancelled 
from the isotropic pressure equation, written in the terms of new variables -u_ and u + Hence, 
we have reduced the second order differential equation to two independent Riccati type 
differential equations, for which exact solutions of the field equations describing the stellar 
structure can be generated by the appropriate choices of the generating function F (R). 



A. Case A: 



F(R) — (3 = constant 

In order to obtain some solutions of the gravitational field equations for our stellar model, 
we assume the generating function F (R) = [3 =constant. In this case equations (|T4*j) take 
the form 

dtl— In dtL-^ In / 1 i-\ 

m_— zm_u_ + u_ = p, —m + — h 2m + u + + u< — p, (15) 

dR dR 

Depending on the sign of the parameter (3, we obtain two distinct classes of solutions of 
Eq. JED 
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(i) A T > 



(ii) A T < 
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where A T = 4 (mi + (3) , e T = =Fl, R± are arbitrary constants of integration. 



B. Case B: 



(16) 



(17) 



F (R) = -2m_« _ + = 2m + « + + v>o + 

Assuming that the solutions of equations (fT4*|) are constants, i. e., u T = uq t =constant, 
we can obtain another solution of the field equations for the static stellar model with a linear 
barotropic equation of state. 

With the use of Eq. (fTTj). equations (fTUj) become 



A 



m + 



dR, 



(18) 



1 
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The explicit forms of the metric functions e A and e v representing the relativistic fluid 
spheres for three distinct classes of solutions are given by 
Case {A) (i) 
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Case (A) (ii) 



e = A n r cos 
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Case (B) 
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where we have denoted a = u^mZ 1 and b = uo + mZ 1 . 

For Case (B), the physical quantities, the energy density p and the fluid pressure p are 
expressed in the following form 
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V = — — — — — — —T^z : — —P = 7P. 



(27) 



8a - 3a 2 - 86 + 2ab - 3b 2 + 2^2 (a - b - 2) (a + b) 
where < 7 < 1. 

Therefore we have obtained an exact solution of the field equations for the perfect fluid 
sphere in isotropic coordinate with a linear barotropic type equation of state. 

By choosing the parameters uo + = 0,Uo- = 2m_ or w _ = 0,«o+ = — 2m + , the metric 
functions given by Eqs. (}2"4*j) and ()25|) generates the line element 



ds 2 = r ±V2 ~dt 2 - r^~ 2 (dr 2 + rW) , (28) 

with the equation of state p = p (stiff or Zeldovich equation of state) for the dense matter. 

It can be shown that the + and - signs are equivalent through the transformation r — ► 
r _1 . Therefore by using the matrix method we have rediscovered this exact solution for the 
stiff fluid matter that was obtained by Haggag and Hajj-Boutros 



To overcome the difficulty of infinite energy density at the center of the sphere, it is 
assumed that the matter distribution has a core of radius ro and constant density po which 
is surrounded by the fluid with pressure equal to energy density. In view of equations (|24|) 
and (|25|) . by choosing b = a we obtain e A ~ e u r^^ b ~~ a \ or equivalently, Uq t = —2a/ 4 =p 2\/2. 
Therefore the interior metric is conformally flat satisfying the zero condition of the Weyl 
tensor. 

a 1 ^ 

By means of the transformation r^2 — ► r, the conformally flat metric becomes 

ds 2 = r^dt 2 2<ir 2 - r 2 dQ 2 . (29) 

(a + V2) 



However, Gurses and Gursey 



In Schwarzschild coordinates Tolman 20] presented this solution over sixty years ago. 



2l\ proved that the Schwarzschild interior metric for a non- 



charged sphere is the unique solution for the Einstein's gravitational field equations which 
is static and conformally flat. 

Surprisingly, the line element (|2U|) contradicts the results suggested by Gurses and Gursey 



IV. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have studied general relativistic stellar models described by 
the line element (jj), and have presented two new classes of solutions, given by Eq. (|20|1- 
()23|) . of the field equations representing the interior structure of the compact objects in the 
framework of general relativity. The class of the solution corresponding to the stiff fluid has 
been discussed in detail. From the physical point of view, the mathematical solutions must 
satisfy certain physical requirements to render them physically meaningful. The following 
conditions or requirements have been generally accepted 

(a) the energy density p and the pressure p should be positive and finite; 

(b) their gradients ^ and should be negative; 

(c) the speed of sound should be less than the speed of light. 

(d) a physically reasonable energy-momentum tensor must obey a trace condition, p > 3p. 

(e) the interior metric should be joined continuously with the exterior Schwarzschild 
metric. 

(f) the pressure p should vanish at the vacuum boundary of the sphere. 
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(g) the structures should be stable under radial pulsations. 

Not many exact solutions for relativistic fluid spheres are known in isotropic coordinates. 
Hence, the main purpose of this paper is to present the matrix method for obtaining the 
exact solutions of the field equations for the compact objects in isotropic coordinates. Eq. 
(|14|) is a first order differential equation in two unknowns w_and u + . To obtain its solutions 
a solution generating function F (R) has been introduced in Eq. (fT3|). Because of the 
mathematical structure of the master equation (|14|). hosts of solutions can be found by the 
appropriate choice of the generating function F (R) or by making an ad hoc assumption for 
one of the variables «_or u + . 

In order to have physically realistic stellar models, one needs to test the solutions obtained 
via this method to ensure that they satisfy all the physical conditions (a)-(g). 
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